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ABSTRACT 

Although the properties of the narrow-line region (NLR) of active galactic nuclei 
(AGN) have been deeply studied by many authors in the past three decades, many 
questions are still open. The main goal of this work is to explore the NLR of 
Seyfert galaxies by collecting a large statistical spectroscopic sample of Seyfert 2 and 
Intermediate-type Seyfert galaxies having a high signal-to-noise ratio in order to take 
advantage of a high number of emission-lines to be accurately measured. 

2153 Seyfert 2 and 521 Intermediate- type Seyfert spectra were selected from Sloan 
Digital Sky Survey - Data Release 7 (SDSS-DR7) with a diagnostic diagram based on 
the oxygen emission-line ratios. All the emission-lines, broad components included, 
were measured by means of a self-developed code, after the subtraction of the stellar 
component. Physical parameters, such as internal reddening, ionization parameter, 
temperature, density, gas and stellar velocity dispersion were determined for each 
object. Furthermore, we estimated mass and radius of the NLR, kinetic energy of the 
ionized gas, and black-hole accretion rate. 

From the emission-line analysis and the estimated physical properties, it appears 
that the NLR is similar in Seyfert 2 and Intermediate-Seyfert galaxies. The only 
differences, lower extinction, gas kinematics in general not dominated by the host 
galaxy gravitational potential and higher percentage of [O iii]A5007 blue asymmetries 
in Intermediate-Seyfert can be ascribed to an effect of inclination of our line of sight 
with respect to the torus axis. 
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1 INTRODUCTION 

Seyfert galaxies are characterized by a bright star-like 
nucleus, whose spectrum shows emission-lines covering a 
wide range of ionization stages. The main components 
of a Seyfert nucleus are: a central source, the broad-line 
region (BLR), a molecular dusty torus, the narrow-line 
Region (NLR) and a possible extended NLR (ENLR). 
The BLR has a typical sub-parsec size and it is sur- 
rounded by an optically thick torus of dust. The NLR 
is characterized by low density clouds {Ne ~ 10^-10® 
cm~^), temperature Tc ~ 10* K, size ~ 1 kpc and full- 
width at half maximum (FWHM) of the emission-lines 
in the range ~ 200-1000 km s"^ The torus confines the 
escaping ionizing radi ation into two o ppositely directed 
cones (see e.g. IStorc hi-Bcrgmann. Mulchaev. fc WilsonI 
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gas is directly ionized by a non-thermal power-law spec- 
trum fr om the central source, ev en if the presence of 
shocks (|Dopita fc Sutherland! !l995! ) and photoionization 
by star- forming regions close to the nucleus (composite 
galajcies, !Hill et abl 1999! *) cannot be excluded. This model, 
called Unified Model and based on the spectropolarimetric 
analys is of NGC 1068, performed by An tonucci fc Mill^ 
(| 19851 ). i s the most accep ted active galactic nuclei (AGN) 
picture (! Antonucci! !l99a ). The Unified Model is able to 
explain the main different characteristics of Seyfert 1 and 
Seyfert 2 galaxies. In the first case, the torus is seen face 
on, then the BLR is directly observable, in the second 
case, the line of sight intercepts the torus and the BLR 
is not visible. Deep CCD images of few nearby galaxies 
made through narrow band interference filters revealed an 
extended struct ure ver y rich of ga s clouds s urrounding the 
activ e nucleus (|Poggel 1988. 1989i: Tadhunter fc Ts vetanovl 
Il989l) . According to ]schmitt fc Kinnev! (|l996! ) Sevfert 1 
galaxies seem to have NLRs with a much smaller size than 
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those of Seyfe r t 2 i f observed along the torus axis. But 
iKraemer et alj (|l998l ) stressed that this cannot be explained 
by a simple orientation effect and t hat NLRs of S eyfert 
1 are physically compact. However, ISchmitt et al.l (|2003| ) 
showed that both Seyfert types have similar distributions 
of NLR sizes, and elongated shapes are observed more 
frequently in Seyfert 2 than in Seyfert 1, whose emission is 
more concentrated toward the nucleus. These results are in 
agreement with the Unified Model if we assume that the 
NLR gas has a disk-like distribution. With the advent of the 
Hubble Space Telescop e fHST), it was possible to observe 
that ionization cones (lAfanasiev et al ] |2007l : IWilson|[l997l . 
and references therein) are in many cases closely asso ciated 
with radio jets and lobes {Wilson fc Tsvetanovlll994l ). and 
dust lanes ^Martini et al.ll2003al lbfr 

If we accept the Unified Model, we should expect that 
the spectroscopic properties of the NLR are similar in 
Seyfert 2, intermediate-type Seyfert and Seyfert 1 galax- 
ies, with an increasing contribution from highly ionized gas 
closer to the central source. Unfortunately, even in case of 
high signal-to-noise spectra with good resolution, the narrow 
components of the recombination emission lines in Seyfert 1 
galaxies are difficult to measure, since they are included in 
and hidden by the broad and bright components emitted by 
the BLR. This is even harder when the application of auto- 
matic fitting procedures is mandatory because of the analy- 
sis of large samples. Spectroscopic differences and similari- 
ties among Seyfert types have been found in visible spectra, 
but o ft en limited to sin gle objects (see e.g. i Kraemer et al.l 
1 19981 : iFraauelli. Storch i-Bcrgmann. & Binettd I2OO0I '). or 
small samples (see e.g. J3enncrt ct al. 2006a, b), while large 
samples are very rarely studied ( Zhang ct a l. 2008). 

The aim of this work is to explore possible differences 
between Seyfert 2 and Intermediate-type Seyfert galaxies, 
through a detailed analysis of a large sample of Seyfert spec- 
tra and taking advantage also of the weak lines, which are 
too often neglected. This work is organized in five sections: 
in Section 2 the selection criteria of the spectra are de- 
scribed; Section 3 describes the measurements of the spectra 
and their classification in Seyfert 2 and Intermediate-type 
Seyfert. The emission-line analysis and the NLR physical 
characteristics are dealt in Section 4 and 5 respectively. The 
last section summarizes the main achievements of this work. 



2 THE SPECTRA SELECTION 

To build a statistical sample of Seyfert galaxies, we 
decided to exploit the Sloan Digital Sky Survey - 
Data Release 7 (SDSS-DR7) database (|Abazaiian et al] 
I2OO9D . We used a diagnostic diagram (hereafter O123) 
based on the [O ll]A3727/[0 lll]A5007 (hereafter O23) and 
[Ol]A6300/[Olll]A5007 (hereafter O13) ratios. The first ra- 
tio gives a measure of the ionization level, while the second 
ratio is a tracer of the hard ness of the ionizing spectru m at 
large radii from the source JShields fc Filippenkol [T990l\ Up 
to now this ratio has never been employed to classify a large 
number of objects. The great advantages of the O123 dia- 
gram are the following: i) the oxygen lines are visible in all 
AGN spectra, ii) oxygen is the only element for which three 
different states of ionization exist in the visible range, iii) 
the employed ratios are sensitive to the ionization parame- 
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Figure 1. Kewley classification l lKewlev et al.|[2006l) on the O123 
diagram applied to the sample of 62036 objects: blue crosses are 
Seyfert 2, red squares are LINERs, green dots are star-forming 
galaxies and yellow triangles are composite galaxies. 



ter and power-law index and are weakly dependent on the 
metallicity and density, iv) they are not contaminated by 
stellar spectral features, therefore it is not necessary to sub- 
tract the stellar template from the original spectra before 
measuring the lines, v) the lines are not contaminated by 
traditional BLR, vi) and finally [Ol]A6300 (and to a lesser 
extent [N ll]A6548,6584) are more sensitive to the hard AGN 
continuum compared with the stellar continuum (of Hll re- 
gions). Nevertheless, this diagram has two main disadvan- 
tages: the ratios are not reddening free and the [O l]A6300 is 
sometimes weak and difficult to be measured especially when 
the signal-to-noise ratio (S/N) of the [Ol]A6300 is less than 
10. From the SDSS-DR7 we extracted all galaxies showing 
[Oii]A3727, [Oiii]A5007, [Oi]A6300 emission- lines, but with 
the constraint S/N([0 1] A6300)> 3, obtaining 119226 tar- 
gets. 

To check the effectiveness of oxygen spectral lines in 
isolating Seyfert galaxies, we carried out a test. First, 
we exploited the database of 85224 galaxies published by 
iKewlev et al.l ((2OO6), finding 62036 galaxies in common with 
our sample. Then, we plotted the Veilleux-Oste rbrock di- 
agnostic diagrams l|Veilleux fc Osterbrock 1987, hereafter 
VO) and we used the relations given in lKewlev et al.l (|2006| ) 
to classify these objects as Seyfert, Low-Ionization Nuclear 
Emission- line Regions (LINERs), star-form ing and compos- 
ite Seyfert-Hll galaxies. Since the lKewlev et al.l (2006) data 
are reddening corrected and the host galaxy spectrum is 
removed, they are not directly comparable with our sam- 
ple. Therefore, we decided to extract from SDSS the spec- 
troscopic information about these already classified objects 
and we plotted the O123 diagram (Fig.[T]). The result was 
extremely encouraging: this plot shows a V-shape where 
the upper side is populated by star-forming galaxies, while 
the lower side is populated by Seyfert galaxies and LIN- 
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Figure 2. Histogram of log(Oi3) values with log{023) in the 
range [—0.1,0]. 



Figure 3. The O123 diagram of our sample of 119226 objects. 
The solid line is the empirical separation described in Eq.[T] 



ERs, which are sufficiently separated. This diagram is not 
new, s ince it was aheady pubhshed bv lShields fc Fihppenkd 
ljl990l ) (see their fig. 13), but the areas indicated by these au- 
thors, where Seyfert galaxies, LINERs and Hll regions are 
expected to be found, do not perfectly match our data. 

From this diagram we derived an empirical separation 
line to isolate Seyfert galaxies. We sampled log(023) in bins 
of 0.1 dex, in the range [—1.2,0.3] and plotting histograms 
of log(Oi3) in bins of 0.1 dex, we estimated the value of 
the minimum between the two peaks of the distributions for 
each bin (as shown in Fig.[2|. Finally, taking into account 
the mean point of the bin, we interpolated these data with 
a polynomial function, obtaining as result the following for- 
mula: 

log(023) = 0.20 - 0.25 X log(Oi3) - 0.39 x [log(Oi3)]^ (1) 

By selecting the objects under the curve we include 97 per 
cent of the original sample of Seyfert galaxies. Of course this 
new sub-sample is contaminated by other objects: in partic- 
ular, it contains 60 per cent of Seyfert galaxies and 31 per 
cent of composite Seyfert-Hll galaxies, but only 4 per cent 
of LINERs and 5 per cent of star-forming galaxies. So, we 
can conclude that about 10 per cent of a sample selected 
by following this method consists of non-Seyfert galaxies. It 
is interesting to note that the lower sequence contains both 
narrow and broad- lined AGN. This is due to the fact that it 
uses only the oxygen lines. In conclusion the O123 diagram is 
useful for the classification of emission-line galaxies in gen- 
eral and not only for narrow emission-line ones. By apply- 
ing Eq.[T]to the whole sample of 119226 objects (Fig.[3]l we 
found that about 16000 populate the Seyfert region on the 
O123 diagram. Since the SDSS spectra are obtained with a 
3 arcsec fibre aperture, a redshift limit 2 0.1 was adopted 
to reduce the flux c ontamination by extra- nuclear sources 
l|Kewlev et al.ll2006h . The lower hmit z > 0.02 is required 
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Figure 4. S/N ratio of the continuum at rest frame 5500 A 
for the 5678 spectra of our sample (left), x'^/dof distribution of 
STARLIGHT fittings to the 5678 spectra (right). 

because the [O ll]A3727 line must be visible in the 3800-9200 
A spectral range covered by the SDSS detectors. 

The total number of galaxies in the Seyfert region that 
satisfy the mentioned conditions is 5678. The left panel of 
Fig. 13] shows the S/N ratio of the spectra measured at rest- 
frame 5500 A. In most cases (91 per cent) it is between 10 
and 40, with a peak around 20. 



3 MEASUREMENTS AND SPECTRA 
CLASSIFICATION 

The VO diagrams are sensitive to the measurement of 
the H/3 emission-line that is made difficult by the pres- 
ence of the underlying stellar absorption. A suitable cor- 
rection is mandatory to avoid an over-estimate of the 
[O lll]A5007/H/3 diagnostic ratio. The fitting and subtrac- 
tion of the underlying stellar continuum from each spectrum 
was carried out with the spectral synthesis code starlight 
(|Cid Fernandes et aLllioOsl . |2007| ) . This code makes a linear 
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combination of synthetic spectra which are reddened and 
then convolved with a broadening function, allowing to ob- 
tain the stellar velocity dispersion. 

Before being analysed with starlight, the spectra 
were processed with iraf. A correction for Galactic ab- 
sorption was first applied to the spectra, by using dered- 
DEN and A(V) extinction values given by NASA/IPAC 
Extragalactic Database (NED). Then, the spectra were 
shifted to rest-frame with newredshift by using z val- 
ues given by SDSS, re-gridded to a dispersion of 1 A/px 
with DISPCOR and converted int o text format. We u sed a s 
base 92 synthetic spectra from lBruzual~fc" Charlod ()2003l l 
by combining 23 ages (from lO'' yr up to 13 x 10^ yr) 
with 4 metallicities (Z = 004,0 008,0.02 and 0.05), and 
ICardelli. Clavton. fc Mathij (jl989l ) (hereafter COM) as ex- 
tinction function. We masked the emission-lines in order 
to improve the quality of the fit. The goodness of the 
STARLIGHT fitting evaluated with is peaked around 1 
(Fig-Hl right panel). The best-fitting synthetic spectrum of 
the continuum for each galaxy was subtracted from the ob- 
served one in order to obtain a pure emission-line spectrum, 
where hydrogen and helium Balmer lines could be correctly 
measured (Fig.[5|. 

The emission-lines were fitted by means of a dedi- 
cated code written in C language. This code, named GGFIT, 
was developed within our group, in order to fit in auto- 
matic way the emission-lines of a group of spectra with a 
multi-Gaussian non-linear least-squares fitting method. In 
the case of the blended lines [Nii]A6548,6584 and Ha or 
[S ll]A6716,6731, the code fits these features simultaneously. 
The output is a table collecting 12 parameters for each spec- 
tral line: 

• the central wavelength (a;o), amplitude {A), sigma (cr), 
and their errors (Aa;o, AA, Acr); 

• the measured fiux (-Fm); 

• the S/N ratio; 

• the parameter, the sum of the square residuals, de- 
fined as = {5x()f + {SAY' + ((5cr)^, where Sxq, 5 A and 5a 
are the differences between the last iteration values and the 
previous ones; 

• the Gaussian fiux (-fc) and its error (AFg); 

• the of the fitting. 

The first step is to determine the line boundaries, i. e. the 
wavelengths at which the line profile fades into the contin- 
uum on both sides of the line. For this purpose, the algo- 
rithm evaluates the intensity of the line profile (/(A)) from 
the peak towards bluer wavelengths until it reaches a mini- 
mum (7min). If Imin < 7c+o"c, where 7c and CTc are the average 
continuum intensity and its error, the associated wavelength 
is assumed as the blue boundary of the line, otherwise the 
algorithm searches for a new minimum. The same procedure 
is repeated on the red side of the line. In order to prevent 
non-convergence, for each line or blend of lines a maximum 
value of the baseline is imposed. Each interval was empiri- 
cally chosen to avoid nearby lines or features and the regions 
where the continuum was evaluated. If this value is reached, 
the algorithm starts again the procedure by comparing /(A) 
with Ic 4- 2ac. If the convergence is not reached, then it tries 
with 3ac. 

Once the boundaries are determined, the flux is ob- 
tained by integrating the spectral feature. The flux error 
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Figure 6. S/N(H^) vs. A(t(H/3) for the measured 5678 spectra. 
The same in the top-right panel, but in logarithmic scale. 



is calculated by means of Ayl and Act. AA is the value of CTc 
determined in proximity of the emission-line, while to esti- 
mate Act a Monte Carlo method is applied during the fitting 
procedure. In particular, starting from the Gaussian param- 
eters, the code generates 3000 possible models near the best 
solution, within the fixed intervals: A — < A < A + CTc, 
CT — 1<CT<CT+1, a;o — l<a;o<a;o + l, and changing ran- 
domly the values of the parameters. A test is performed 
and only the models with a significance level above 70 per 
cent are considered. This procedure was applied only to fea- 
tures with one or two lines, for computational time reasons. 
Ayl is always lower than CTc, whereas Act and Aa;o both de- 
pend on the S/N. In Fig.[6]the case of H/3 is presented. As it 
can be seen. Act oc (S/N)". Following Corsini ct al. (199^) 
and by applying a linear fitting to log Act vs. log S/N, a sim- 
ple relation to estimate Act can be derived. In this case. Act 
is given by: 



log(ACT) = -1.21og(S/N) +0.7 



(2) 



but similar results are obtained for all the measured lines. 
Aio is similar to Act and it correlates well with Act and 
S/N. The correlation coefficient of the interpolation is 0.95. 
The typical wavelength error is about 0.5 A when the S/N 
is around 10. 

In order to check the reliability of the automatic fit- 
ting procedure, about 1300 spectra were measured by hand. 
The line boundaries limits and the fiuxes obtained are the 
same with very rare exceptions. The number of failures in- 
creases when the S/N < 5, therefore only the lines with a 
S/N > 5 were measured. They are listed in Table[T] with the 
maximum number of components adopted in the fitting, the 
ioniz ation potentials and the logarithm of the critical den- 
sity (|De Robertis fc Osterbrocki Il986l : iMendozal Il983l '). No 



other line was taken into consideration because it was too 
weak to be detected, as [Fex]A6374, He and higher order 
Balmer lines, or in blend but presumably giving a small con- 
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Figure 5. An example of correction for stellar component. The original spectrum is plotted in black, the fit to the stellar component in 
red and the residual pure emission-line spectrum in black, dotted line. 
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Figure 7. Examples of lines fitting. The observed spectrum 
(black dots) is overlapped to the results of the fitting (solid line). 
Residuals are plotted on the top of the diagrams (dotted line). 



tribution to the total flux of the blended line, as [S lll]A6312, 
[Ca v]A6087 (|Nagao et al.ll2003l ). Our software was able to fit 
successfully 4254 spectra: the others were rejected because 
either the most important emission-lines (i.e. those used in 
diagnostic diagrams) had S/N < 5 or (in a few cases) their 
fitting did not converge. Fig.[7] shows an example of the fit- 
ting. 



The classification of Seyfert galaxies is based on 
the detection of the Ba lmer line broad components 
(lOsterbrock fc FerlandllioO^ ). but this criterion is subjective 
and it depends on the spectral quality. Additionally, visu- 
ally classifying a large number of spectra is very difficult 
and extremely time-consuming. In this work the criterion 
adopted for the spectral classification is based on the pres- 
ence of a second component in the Ha emission-line and on 
its FWHM. Gclbord ct al. (2009) claimed that a good fitting 
of Hq should be obtained using three components, a narrow, 
an intermediate and a broad one. However, it is in general 
difficult to find the corresponding components in H/3, which 
is fundamental in our analysis. Therefore, in this work we 
decided to fit both H/3 and Ha with two components, one 
narrow and one broad. Out of the 4254 measured spectra, 
1938 show a secondary Ha component. In Fig.[5]we present 
the FWHM(Ha) distribution of both components after cor- 
rection for instrumental width (R=1800). The FWHM(Ha) 
distribution of the second component shows two peaks, one 
between 500 and 1000 km s~^ and the other between 1500 
and 2200 km s"\ respectively. If the FWHM(Ha) of the sec- 
ond component is lower than 2000 km s~^ the line is hidde n 
in the Ha+[N ll] A6548,6584 feature dCelbord et al] bOOgl V 
In this case, the SDSS spectral resolution is not sufficient to 
separate the true components. Therefore, since we could not 
be sure whether the second Ha component is a real broad 
line or a mathematical result of the fit, we decided to ex- 
clude the broad components with FWHM < 2000 km s~^. 
In addition, we compared the widths of Ha and H/3 lines 
for all 4254 measured spectra (Fig. [9]). We noticed that the 
narrow components of both lines are well distributed around 
the 1:1 line, while the broad components are in good agree- 
ment only when the FWHM(Ha) > 2000 km s'^ Below 
this value, the widths of the two emission-lines are totally 
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Table 1. Measured lines. 



line 


comp 


IP low cV 


IF nign eV 


log(Afc cm ^) 


[O ii]A3727 


1 


13.62 


35.12 


3.7 


[Neiii]A3869 


1 


40.96 


63.45 


7.0 


[Neiii]A3967 


1 


40.96 


63.45 


7.0 


[S ii] A4070 


1 


10.36 


23.34 


6.4 


H(S A4102 


1 


13.60 






H7 A4340 


1 


13.60 






[0 ni]A4363 


1 


35.12 


54.94 


7.5 


He II A4686 


1 


24.59 


54.42 




[Ar iv]A4711 


1 


40.74 


59.81 


4.4 


[Ariv]A4740 


1 


40.74 


59.81 


5.6 


H/3 A4861 


2 


13.60 






[O iii]4959 


2 


35.12 


54.94 


5.8 


[O iii]A5007 


2 


35.12 


54.94 


5.8 


[Ni]A5199 


1 





14.53 


3.3 


[Fevii]A5721 


1 


99.1 


124.98 


7.6 


HeiA5876 


1 





24.59 




[Fevii]A6087 


1 


99.1 


124.98 


7.6 


[O i]A6300 


1 





13.62 


6.3 


[O i]A6363 


1 





13.62 


6.3 


[Nii]A6548 


1 


14.53 


29.60 


4.9 


Hq A6563 


2 


13.60 






[Nil] 6583 


1 


14.53 


29.60 


4.9 


[S ii]A6716 


1 


10.36 


23.34 


3.2 


[S ii]A6731 


1 


10.36 


23.34 


3.6 


[Ar iii]A7136 


1 


27.63 


40.74 


6.7 


[O ii]A7320 


1 


13.62 


35.12 


6.8 


[O ii]A7330 


1 


13.62 


35.12 


6.8 



Notes. The table shows the measured lines, the number of com- 
ponents used in the fitting, their ionization potentials (IP low) 
and the potential of the following ionic species (IP high), and the 
logarithm of the critical density. 




log FWHM (Ha) 

Figure 8. Distribution of FWHM(Ha) for the first compo- 
nent (solid line) and the second component (dashed line) of the 
emission-line. 
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Figure 9. FWHM(H/3) vs FWHM(Ha) for the 4254 measured 
spectra. Crosses are the first component, squares the second com- 
ponent. 

diflterent. In particular, the second component of H/3 could 
not be a real broad line, or it could be poorly fitted, be- 
cause too weak. When FWHM(Hq) < 1250 km s"^ the sec- 
ond component of H/3 disappears. Therefore, all the spectra 
with narrow Ha or with a second component narrower than 
1250 km s~^ were classified as narrow-line spectra. The er- 
rors were taken into account, and the condition was fixed to 
FWHM(Ha) - 3 x AFWHM(Hq) < 1250kms~\ The Ha 
flux was obtained as the sum of two fiuxes when a second 
component was present. The spectra showing a second Ha 
component with FWHM(Ha) -f 3 x AFWHM(Hq) > 2000 
km s~^ were classified as broad line spectra. In this sec- 
ond class we imposed an additional constraint: the reliabil- 
ity of the broad component was confirmed only when the 
height of its profile above the continuum at rest-frame 6531 
and 6595 A was larger than 2a^, where ctc is the root mean 
square of the continuum intensity. The factor 2 was adopted 
based on visual tests on the detection of the broad compo- 
nent. Fig. [To] shows an example of broad- line spectra which 
satisfy the previously mentioned conditions. The distribu- 
tions of the second component FWHM(Hq) of these groups 
are showed in Fig. 1111 Finally, 598 objects were classified as 
broad emission-line galaxies and 2836 as narrow emission- 
line galaxies. Then, we applied t he VQ diagno s tic di agrams 
and the relations introduced bv lKewlev et al.l (|2006l ). Only 
the objects belonging to the Seyfert region within la error 
in the diagnostic ratios were considered. We obtained two fi- 
nal samples of 521 Intermediate-type Seyfert (S-I) and 2153 
Seyfert 2 (S2) galaxies. 

As a final check w e plotted the l o g([0 i ii]A5007/Hq;) 
vs log([Nii]A6584/Ha). ICelbord et all lj2009h found that 
this diagram is able to separate the S2 from the 
other Seyfert galaxies. They found that S2 galaxies have 
log([Oiii]A5007/HQ) > and log([N ii]A6584/Hq) > -0.5, 
while our boundaries are log([0 ill] A5007/Hq) > —0.6 
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Figure 10. Example of H/9 (top) and H« (bottom) with a broad 
second component. Tlie observed spectrum is in black dots, the 
components are in dashed line and the residuals in dotted line. 
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Figure 11. FWHM(Ho) distribution for the second component. 
The spectra showing Ha second component with FWHM < 
1250 kms~^ are in solid line, while those having FWHM > 
2000 kms~^ are in dashed line. 



and log([Nll]A6584/HQ) > -0.3 (Fig.[T21). On the other 
hand, the correspondence is very good for the S-I galax- 
ies. This discrepan cy could be due to the fact that 
iGelbord et al.l l|2009l ) considered only spectra showing for- 
bidden high-ionization line (FHIL), and so it is likely that 
[O lll]A5007/Ha ratio is higher than the mean value we found 
in S2. 
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Figure 12. The logarithm of [N ll] A6584/H« ratio plotted against 
the logarithm of [O ill] A5007/Ha ratio. S2 galaxies are indicated 
with red plus and S-I galaxies with blue squares. 
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4 EMISSION-LINE ANALYSIS 

In the following sections, we present the statistical analysis 
of the emission-line ratios. We looked for the correlations by 
means of the Spearman rank correlation coefficient (rs) that 
is in general adopted when the distributions are a priori 
unknown. We used the two-tail test to verify the null hy- 
pothesis of each correlation, that is the probability that the 
correlation does not exist, and we assumed q = 0.1 per cent, 
as significance level. We calculate d the t paramet er of the 
Student distribution by following (|Wall et al.ll2003l ). and we 
compared those values with tabulated critical t- values {tc-it) 



Figure 13. A(y) distribution for the samples: S2 (red solid line), 
and S-I (blue dashed line). 

for a given number of degrees of freedom and a significance 
value. When t > tcrit, the correlation is significant at the 
chosen level and the null hypothesis can be rejected. 

To compare samples, we used the Kolmogorov-Smirnov 
(K-S) test, to verify when two samples are drawn by the 
same population. When the observed parameter A, that 
is the maximum difference between two distributions, is 
greater than the theoretical value D, the two samples are 
different with a significance level of 1 per cent. 
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Table 2. Correlation coefficients between lines of similar ioniza- 
tion potential. 



Figure 14. Stellar reddening (A^^"^) versus gas reddening 
Symbols are like in Fig. 1121 



4.1 Internal reddening 

The internal reddening correction was calculated using the 
parametrization introduced by CCM under the assumption 
Hq/H^ = 3.1 (|Osterbrock fc Ferlandl [20061 ). Fig.[T3] shows 
the distributions of the derived values of A(l/). The median 
values are 1.4 (90 per cent of the sample has A(V) < 2.5) for 
the S2 and 0.9 (90 per cent of the sample has K{V) < 1.8) 
for the S-I. Assuming a level of significance of 1 per cent, 
the Kolmogorov-Smirnov (KS) test gives 0.3, clearly dif- 
ferent from the theoretical value of 0.08, indicating a sta- 
tistical difference of A(V^) distribution between S2 and S- 
I. This result is expected in S-I and Seyfert 1 galaxies, 
because we are looking deeper inside the NLR, and the 
dust quantity is p resuma bly lower because of the winds 
l|Calzetti. Kinnev~fc Storchi-BerB ; mannlll996t ) and the dust 



sublimation ( Mullanev et al. 20091) 



The A(l/) values appear weakly correlated with stellar 
absorption (A(K)), evaluated by starlight (Fig.[T3}. In- 
deed the Spearman rank correlation coefficient is r^ = 0.55, 
but the correlation is significant since t = 33.3 >> tcrit = 
3.3 for N = 2557 targets (S2-f S-I). 

Our dis t ributi on is very similar to that obtained 
by IGu et all (|2006l ). who measured the extinction of 65 
nearby Seyfert 2 nuclei from II7/II/3 ratio. A^iVi,) val- 
ues are distributed in a small interval and are very 
similar for both the samples. A large difference be- 
twe en the stellar and gas extinction was a lso found 
by ICalzetti. Kinney fc Storchi-Bergmarml (| 19941 ) and by 
ICalzettil in starburs t galax ies. The same result was 

obtained by Bennert et al.l (|2006al ) in a detailed analysis of 
the NLR in NGC 1386. This would suggest that most of the 
dust is mixed with the gas. 



A5S76 1 
A4886 1 
X3S69 1 
A7136 1 
X3727 1 
X5199 1 
A6300 1 
A65S4 ' 



A46S6 
A5007 
X6007 
A5007 
A6724 
A6724 
A6724 
A6724 



154 

354 
1126 

225 
2110 

165 
2105 
2100 



0.17 
0.60 
0.77 
0.75 
0.58 
0.38 
0.76 
0.51 



2.1* 
14.1 
40.5 
16.9 
32.7 
5.2 
53.6 
27.2 



133 
199 
431 
114 
513 
88 
508 
509 



0.38 
0.49 
0.77 
0.78 
0.71 
0.50 
0.80 
0.69 



4.7 
7.9 
25.0 
13.2 
22.8 
5.4 
30.0 
21.5 



Note. The asterisk indicates that the correlation is not signifi- 
cant. 



4.2 Line intensity correlations 

The emission-line fluxes normalized to H/3 a n d cor rected 
for reddening were analysed following iKoskil (| 19781 ). The 
emission-line ratios were grouped according to high and 
low ionization lines (Table [2| and plotted vers us the 
[OllllA5007/H^ or [S ll]A6724/H^, respectively. As in lKoskil 
(1978) we added the He: A5876/H/3 vs. Hen A4686/H^ di- 
agram, in order to compare the helium recombination lines. 
These diagrams show the correlations between lines with 
similar ionization states (Figures ll5lll6|) . Although the mean 
S/N of our spectra is lower than that of Koski, we found 
similar results by plotting a larger sample of objects. 

Table[2]shows the number of targets (N), the Spearman 
rank correlation coefficient (ts) and the parameter of the 
two-tail test (t). The critical value for the minimum number 
of targets, N = 88, is tcrit = 3.4, therefore all correlations are 
significant at 0.1 per cent, but the first one for S2 galaxies, 
indicated by an asterisk. 

[Nelll]A3869 and [Arlll]A7136 lines show a good corre- 
lation with the [O in] A5007 fine in both samples. He 11 A4686 
is well correlated only in S2 galaxies, while in S-I galaxies the 
correlation is significant, but weaker, likely because of the 
presence of the underlying broad component. Good or fairly 
good correlations exist between [Oi]A6300 and [Nii]A6584 
vs. [Sll]A6724, while [Nl]A5199 vs. [Sll]A6724 is weak in 
both samples, probably because [Nl]A5199 is very close to 
the Mg I A5175 triplet and its measure can be strongly af- 
fected by the stellar continuum subtraction. We found also 
good correlations between lines with similar critical density: 
[Sll]A6724 vs. [Oll]A3727, [Arlll]A7136 vs. [Oll]A7325 (only 
for S2) and [Ne ill] A3869. 

On the other hand, there is no correlation between the 
He I and Hell lines, and [Arlll]A7136 vs. [Ol]A6300 and 
[Sll]A4074. These last lines have similar critical densities, 
but very different ionization potentials, supporting the hy- 
pothesis of two independent ionization zones. We do not find 
any correlation between the [Olll]A5007 line and the high 
ionization fines [Fe vil]A5721,6087 and [Ar lv]A4711,4740. 
The [Fevil] lines have an higher critical density compared 
to the [Olll] lines, while the [Ariv] lines have higher ion- 
ization potential and lower critical density. Therefore, these 
lines could be formed in different regions compared with 
the [Olll]A5007 line. [Fe vil]A5721,6087 were detected in few 
spectra, but the percentage increases from S2 to S-I (see 
Table [Sj . These lines are formed in the coronal-line region 
(CLR) an d are very important in t he an alysis of the NLR 
structure. iMuravama fc Taniguchil (|l998l ) claimed that in 
the context of the Unified Model three kinds of CLR exist: 
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Figure 15. Koski diagrams for high ionization hnes (see text). The fluxes are relative to H/3 and reddening corrected. Symbols are like 
in Fig.[T2l 
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Figure 16. Koski diagrams for low ionization lines (see text). The fluxes are relative to H/3 and reddening corrected. Symbols are like 
in Fig.[T2l 



those associated with the torus, with the NLR and with the 
ENLR. The condition to form the coronal li nes is a dust free 
medi um with densities of 10^-10*'^ cm~^ () Ferguson et al.l 
1 19971 ). If the CLR is far away from the source, it is likely 
to ha ve a very low density, ~ 1 cm~^ (|Korista fc Ferlandl 
1 19891 ). 



Since the critical density of the [Fevn]A6087 line is 
very similar to that of the [O l ll]A43 63 line, we followed 
iNagao. Muravama. &: Taniguchil (|200ll ) and analyzed a pos- 
sible correlation between these two lines. About 50 per 
cent of the spectra where the [Oni]A4363 line is detected 
shows the [Fevn]A6087 line. In particular, we calculated 
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Table 3. Number of spectra with [Fevil] and [Olll]A4363 lines 
and their percentage with respect to the samples. 



samples [FcVll]60S6 [O III] A4363 [O ill] A4363+ [Fc Vll] A6087 

S2 96 (4 per cent) 86 (4 per eent) 38 (2 per cent) 

S-I 97 (19 per cent) 186 (36 per cent) 73 (14 per cent) 
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Figure 17. The logarithm of the [Fo vn]A6087/[O ni]A5007 ratio 
is plotted against the logarithm of the [O III] A4363/[0 III] A5007 
ratio for the two samples. Symbols are like in Fig. 1121 

[Fevii]A6087/[Oiii]A5007 and [O in]A4363/[0 ni]A5007 ra- 
tios. We did not find correlation for S2, r3=0.28 {t — 1.8 < 
tcrit ~ 3.6, N=38) and a weak but significant correlation 
for S-I, r,=0.47 {t = 4.5 > tc„t = 3.4, N=73) (Fig.[T3. 
The low value of the correlation index in the S-I sam- 
ple could be due to larger errors in the measurements of 
the flux of [Oin]A4363, because of the presence of a non- 
negligible broad H7 component. It is interesting to note 
that the distribution extends to higher values in the S-I 
case. Indeed, if the temperature increases and the density 
decreases, then the [O iii]A4 363/[0 ni]A5007 ratio increases 
(lOsterbrock fc Ferlandll2006l ). This is consistent with the as- 
sumption that the coronal lines are forming in a medium 
with low or intermediate density and high temperature. The 
same critical density is not a sufficient condition to support 
the idea that these lines are emitted in the same region. 
Their ionization potentials are very different, therefore it 
is likely that these lines are in general emitted in different 
regions and in some cases in partially overlapped regions. 



5 PHYSICAL CHARACTERISTICS 
5.1 Densities and Temperatures 

The values of density and temperature were calculated by 
means of the IRAF task temden. The line ratios used 
to determine the densities were [S n]A6716/6731 (hereafter 
RS2) and [Arlv]A471 1/4740 (hereafter RAr4), while the ra- 
tios used for the temperature were [O in] (4959-f 5007) /4363 
(hereafter R03), [S n](6716+6731)/(4068-f 4076) (hereafter 



Table 4. Median values of electron temperature for the samples. 
The meaning of D and A is explained at the beginning of Section 
4. 



temperature 


S2 


S-I 


N S2 


N S-I 


D 


A 


T02 


12000 


10000 


68 


13 


0.49 


0.14 


TS2 


11000 


15000 


63 


53 


0.30 


0.30 


T03 


17500 


25000 


84 


170 


0.21 


0.42 



RS2t), [O ii]A3727/A7325 (hereafter R02). A temperature of 
10000 K was assumed for the density determination. Fig. 1181 
shows the distribution of the densities calculated by means 
of RS2. This ratio was measurable in all spectra, but rea- 
sonable values (RS2 > 0.29, theoretical lower limit) and so 
reliable estimates of the density were possible in 2316 out 
of the 2674 objects, 1846 S2 and 470 S-I (87 per cent of the 
total sample). The median value is ~ 250 cm"'', most of 
objects (~ 2300 galaxies) show values lower than 500 cm"'', 
while only in 97 cases we observe densities higher than 1000 
cm""'. Only in 11 objects it was possible to estimate the den- 
sity with argon lines, 7 S2 and 4 S-I. The obtained values are 
of the order of 10^ cm~^ (Fig. [19}. This is not unexpected, 
since the [Ar iv] transitions are characterized by higher crit- 
ical density values than those of the [S 11] lines, and therefore 
they are probably emitted by gas in different physical condi- 
tions. Interestingly, the [Olll]A4363 transition, which shares 
a similar ionization potential, has an even higher critical 
density, and this means that the gas emitting this line is 
also characterized by a larger electron density, likely around 
10'* cm"'' or even more. Moreover, both lines [Arlv]A4711 
and [Arlv]A4740 are generally weak, also after reddening 
correction, showing median intensities normalized to H/3 of 
0.07 and 0.09, respectively. Since most of the spectra, where 
they could be detected and measured, have S/N ratio be- 
tween 20 and 40, we can in principle conclude that these 
lines are visible only in few galaxies, because spectra with 
high S/N ratio in the continuum are mandatory. Anyway, 
this could be only a necessary but not sufficient condition. 
Indeed, 66 per cent of our sample has S/N > 20. More- 
over, one can speculate for example that [Oni]A4363 has a 
median intensity only twice higher than [Arlv]A4711 and 
[Arlv]A4740, and it was detected in 272 spectra having S/N 
ratio between 15 and 40, notwithstanding the fact that oxy- 
gen is about 200 times more abundant than argon, assuming 
solar abundances. Therefore, it is also possible that the gas 
of the NLR is made essentially by a low and a high density 
medium, and that the second one has the necessary electron 
density to coUisionally populate the oxygen auroral line and 
simultaneously suppress the [Ariv] lines. 

The estimate of temperature can be difficult because 
sometimes the lines involved are weak or even very weak. 
Anyway, reliable estimates were obtained with R03 in 254 
objects, 84 S2 and 170 S-I, with R02 in 81 objects, 68 S2, 
13 S-I, and finally with RS2t in 116 reliable measures, 63 
S2 and 53 S-I. Fig. 1201 shows the distributions of the so de- 
termined temperatures for both samples while the median 
values are reported in Table [J] We used the electron density 
obtained with RS2 as input to temden. Through a two- 
sample Kolmogorov-Smirnov test (see Table |4)| it appears 
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Figure 18. Density distributions calculated from RS2 (see text). 
Colours and lines are like in Fig. 1131 




Figure 19. Density distributions calculated from RS2 (thick line) 
and RAr4 (dashed line). 



that the distributions of the temperatures derived from [O ll] 
and [S ll] are very similar for S2. On the other hand, the tem- 
peratures obtained from [O III] have a completely different 
distribution. In principle, this result could be indicative of 
the presence of ion stratification according to the ioniza- 
tion potentials. In contrast, it could be simply the effect of 
the use of A'e([Sn]), which indicate too low electron den- 
sity values to coUisionally pump the [O in]A4363. The use of 
A'e([Arlv]) should be more correct, nevertheless the critical 
densities of [Ariv] lines are lower than that of [Olll] A4363, 
and this suggests that A'^o([Ar iv])=10^-10'* cm"^ could be 
a lower limit for the density of the gas emitting the [O ill] 
auroral line. 

By assuming that photoionization from the active nu- 
cleus is the main physical mechanism taking place in the 
NLR, we can conclude that [O ill] A4959, 5007 are likely emit- 
ted both by a low density medium at higher temperature 
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Figure 20. Distribution of temperature in logarithmic scale, 
measured by means of different line ratios: Tc([Olll]) (top), 
rc([Oll]) (middle), Te([Sll]) (bottom). Colours and lines are like 
in Fig.[T51 



reflecting ion stratification, and by high density clouds with 
high ionization, where the physical conditions are such as 
to allow the [O lll]A4363 to form and become detectable and 
measurable. The high values of temperatures measured with 
[O III] lines in the S-I sample, seems to indicate that we are 
observing the part of the NLR closer to the AGN. However, 
we stress that this result could be an effect of non perfect 
subtraction of H7 broad component. 

5.2 Ionization parameter 

The ionization level is defined by the ionization parameter 
U, which is the ratio between the flux of ionizing photons 
and the hydrogen density, 

U = (3) 

where Qion is the number of ionizing photons per sec- 
ond emitted by the source. The ionization parameter 
could be eval uated u s ing t wo ratios, [O iii]A5007/H/3, as 
sugg ested by ICohenI (1 19831). and the most recent ratio 
[O ii]A3727/[0 iii]A5007 iKomossa fc SciiuSl 19971 ). We plot- 
ted the [Oll]A3727/[Olll]A5007 vs [O lll]A5007/H^ diagram 
(Fig. l21|l . The trend is clear in both the samples: r^ = 
-0.65 {t = 39.7, N=2153) for S2, r^ = -0.56 (t = 15.4, 
N=521) for S-I. Cohen & Osterbrock (1981) defined high 
ionization when [O ill] A5007/H/3 > 10. They analyzed dif- 
ferent line ratios in order to define a non-arbitrary way 
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Table 5. Two-sample Kolmogorov— Smirnov test concerning the 
ionization potential distributions with 1 per cent of significance. 



samples 


D 


A 


S2 ^ S-I 


0.08 


0.38 


(S2 ^ S-I)a4686 


0.14 


0.25 


(S2 ^ S-I)^6087 


0.23 


0.21 



Notes. (S2 -o- S-I)a4686 ^iid (S2 S-I)a6087 indicates the com- 
parison between the samples when the spectra show respectively 
the Hell A4686 and [Fevil]A6087 emission-lines. 



to classify high and low ionization level. In their fig. 2, 
[Fevil]A6087/H/3 vs. [O lll]A5007/H/3 diagram, they found 
that [Fe vil]A6087/H/3 shows a cutoff: [Fe vil]A6087 is not ob- 
served, when [Olll]A5007/H^ < 10. From Fig.[2T]we can es- 
timate that this limit corresponds to [O ii] A3727/[0 ill] A5007 
lower than 0.5 and 0.4 respectively for S2 and S-I. It cor- 
responds to a log(U) greater than -2.5 and -2.4 respec- 
tively for S2 and S-I, by assuming the relation between the 
[O n]A3727/[0 ll l]A5007 ratio a n d the ionization parameter 
U introduced bv lPenston et all (|l990h : 



log U = -2.74 - log(J3727//5007). 



(4) 



Using Eq.[3]we estimated the ionization parameter for both 
the samples, and we applied the KS test subdividing each 
sample in two sub-samples of spectra (Table [5]), the first 
one showing the Hell A4686 emission-line and the second 
one showing the [Fevil]A6087 emission line. There is a sig- 
nificant difference between S2 and S-I galaxies when all the 
spectra are taken into account. It is interesting to point out 
that this difference becomes less significant when we com- 
pare the objects emitting Hen A4686 and disappears when 
the [Fevil]A6087 is taken into account (Fig.[22Il- The ioniza- 
tion parameter depends on the luminosity and is inversely 
proportional to density and distance from the source. The 
different distributions suggest that one or more of these pa- 
rameters change. If the luminosity changes, then S-I and S2 
are intrinsically different objects and the Unified Model fails. 
On the other hand, since the average electron density seems 
to be the same in the two samples (see Section [5]), changes 
in the ionization parameter could be a distance effect. The 
ionization is highly stratified in the NLR and in S-I galaxies 
we are able to observe closer to the active nucleus, therefore 
we can speculate that S2 galaxies showing Hen and [Fevil] 
have a wider aperture angle of the torus, which is thinner in 
the direction orthogonal to the line of sight. Indeed, by plot- 
ting log U vs. A(V^) we can observe that S2 galaxies showing 
[Fevil] are characterized by a relatively low extinction, with 
values similar to those observerd in S-I galaxies fFig. [23)) . 



5.3 Ionized gas mass and radius 

[h3 l|2009l ) suggests that the mass budget of the NLR can 
be accounted by mass loss from evolved stars. From Ha 
(or H/3) l uminosity it is possible to d etermine the ionized 
gas mass jOsterbrock fc Ferlandll2006l ). The author found 
Mnlr ~ 3 X 10" M0 in a region of 200 x 400 pc. This value 
comes from the population statistics of the Palom ar sur- 
vey assembled in IHo. Filippenko. fc Sarge"n3 (|2003l ) and it 



g 0.5 - 



0.5 - 




0,5 1 
log (|0 111|A5007/H^') 



L5 



Figure 21. The logarithm of the [O III] A5007/H/3 ratio is plotted 
against the logarithm of the [O ll] A3727/[0 III] A5007 ratio. Sym- 
bols are like in Fig. 1121 




Figure 22. Ionization parameter distributions. Top: S2 sample is 
indicated in red solid line and the S-I sample in blue dashed line. 
Middle: the distributions taking into account only the spectra 
with the Hen A4686 measured, bottom: the distributions taking 
into account only the spectra with the [Fevil]A6087 measured. 



probably underestimates the mass by no more than a fac- 
tor of 2. Comparing this mass with the amount of mate- 
rial shed by evolved s t ars, d etermined by the models of 
iPadovani fc Matteuccil (119931 ) . Ho found ~ 0.05 Mq 
yr~^. Therefore, stellar mass loss can sustain the gas reser- 
voir in the NLR if the stellar debris survives for a lapse 
of time longer than ~ 2 x 10^ — 10^ yr before it dis- 
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Table 6. NLR mass (Eq.jSJl and radius (Eq.[7J distributions for 
three different A^c- 



-1 



Figure 23. The logarithm of the ionization parameter is plotted 
against the reddening derived from emission-lines. Symbols are 
like in Fig. ll3l Black filled circles are Seyfert 2 showing [Fevil], 



sipates and merges with the surrounding hot interstellar 
medium. The mass loss from evolved stars accounts for the 
high metallicity found in the ionized gas of Seyfert galax- 
ies (jStorchi-Bergmann et aP 1 19981 : iHamann et aD l2002h . 
Indeed, only in rare cases the metallicity appears 



to be sub-sola r llGroves. Heckman. fc KauffmannI 20061: 
llzotoy fc ThuanI l2008l l. Following 'Osterbrock fc Ferland' 
l|2006l ), we determined the ionized gas mass by means of 
the Eq.El 



M = 6.67 X 10 



^33 Lhp 
iVc 



Mr: 



(5) 



By assuming the densities derived from the RS2 ratio, we ob- 
tained a mass distribution having a median value of 10® Mq 
and with 80 per cent of the sample having a mass lower than 
3.4 X 10® Mq. Deriving the NLR radius under the assump- 
tion of spherical distribution, we found a median value of 20 
pc, with 80 per cent of the sample under 40 pc. Of course 
these values are too low compared with the NLR radii. Even 
if a filling factor of 10"'^ is assumed, the median radius is 200 
pc. The 80 per cent of the sample has a radius lower than 400 
pc, a more reasonable value but still rather low compared 
with estimated values typically found in the literature (e.i 
iBennert et al. I l2006bl : iKraemer et"al]|2009l : lYoneharall2006 
These results suggests that the NLR mass calculated with 
RS2 density values is probably underestimated. If we as- 
sume a low density medium (e.g. 1-10 cm~'^), the median 
value of the mass becomes about 2.5 x 10^~* M©, and the 
median values of the radius changes between 200 and 900 
pc, more in agreement with the observations. This simple 
exercise supports the idea that most of the NLR is made of 
a low density medium. Alternatively, it is possible to deter- 
mine the gas mass from the me asure of the reddenin g and 
from a given dust-to-gas ratio (|Fu fc StocktonI l2007l ). As- 
suming the Galactic dust-to-gas ratio we obtain the Eq.[6l 



No cm -3 


M (Mq) 


mean (pc) 


median (pc) 


<80 per cent (pc) 


S2 


RS2 


lo'^ 




170 


140 


240 


10 


2.5 


10^ 


750 


680 


950 


1 


2.5 


10^ 


2400 


2200 


3000 


S-I 


RS2 


lo'^ 




200 


150 


260 


10 


2.5 


10^ 


1000 


850 


1200 


1 


2.5 


10^ 


3100 


2700 


4000 



M = 5.28 X 10^^ cm"^ mag"^ E(B ~V)mpQ d\ (6) 

where Q. is the solid angle subtended by the NLR and dA the 
angular distance. Unfortunately it is not straightforward to 
apply this formula in our case. The fiber diameter is obvi- 
ously fixed and it is not possible to get Q. from observational 
measurements. In principle one can assume fixed dimensions 
of the NLR, but this is not a good idea because the typical 
NLR radius varies in a wide range (100 to 1000 pc or more). 
Nevertheless if we determine the NLR mass from Eq.O Eq.[6] 
is useful to estimate the NLR dimensions. In fact by com- 
bining these equations we can express Q,d\ as & function 
of measurable quantities. This gives us the projected NLR 
surface and then the NLR radius 



7? = 0.16 



M 
k(V) 



pc 



(7) 



where tha mass M is given in solar masses. Table [6] contains 
the radius distributions, for both the samples, obtained with 
the three calculated mass median values through Eq.[7l In 
order to obtain the same radius distribution with both the 
approaches (Eq.[S]and[7]) and with 10® M© we have to impose 
a filling factor of ~ 0.001 — 0.01 when the classical approach 
is used. 

As already pointed out by IFu fc Stockton! (|2006| ). and 
in agreement with what we said above, we can reasonably 
assume that most of the emitting volume of the NLR is dom- 
inated by a massive and highly ionized diffuse medium at 
very low density (1-10 cm~'^), containing less ionized clouds 
with density of ~100-500 cm~^, and likely clouds and/or 
filaments at higher density (~ 10^-10^ cm~'^). 

In order to compare the NLR radii and masses of the 
two samples, we assumed a density of 10 cm~'^. The NLR 
radii fFig. l24p are slightly larger in the S-I sample, likely 
because of lower A(V) values measured in these objects 
(Fig. ll3|) . On the other hand, the mass distributions are simi- 
lar for both the samples (Fig. l25p . The median value is about 
3 X 10^ M© , and 90 per cent of the sample has a mass lower 
than 10* Mp). 



5.4 Kinematics 

The measured FWHM (FWHMm) was corrected for instru- 
mental broadening to obtain the intrinsic value (FWHMi). 
As instrumental widt h we used th e nominal SDSS spectral 
resolution, R = 1800 (|Burles et al., 1999) which corresponds 
to 167 km s~^. 



FWHh'h = ^/{FWHMm X c/Ai)2 - 167^ kms" 



(8) 
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NLR radius (pc) 



Figure 24. NLR radii distributions using A^c = 10 cm ^ . Colours 
and lines are like in Fig. ll3l 



Figure 26. The distributions of the mean FWHM (km s-i) of 
the emission-lines measured in each target. Lines and colours are 
like in Fig.[T3l 
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Table 7. Correlation coefficients between the logarithms of 
FWHMetar and FWHM„as 
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Figure 25. Ionized gas mass obtained by Eq.[5] with Ne 
10 cm^ . Lines and colours are like in Fig.lTsl 



Since a low S/N ratio makes the line profile difficult 
to be measured, we considered only tlie briglitest lines, 
namely: H/?, Ha, [Ol]A6300, [Nll]A6584, [S ll] A6716,6731, 
[Neiii]A3869, [Oiii]A5007 and Hen A4686. Furtfiermore, we 
calculated the mean FWHM for each object (see Fig. l26p . 
using the previously mentioned lines. A very little differ- 
ence between the samples is found by means of the KS 
test (D=0.08, A = 0.15, significance=l per cent). The me- 
dian values are about 270 km s^^ and 300 km s^^ respec- 
tively for S2 and S-I, with the 80 per cent of the sam- 
ples showing FWHM lower than 360 and 390 km s~^ re- 
spectively for S2 and S-I. Concerning a possible connec- 
tion between gas velocity dispersion and stellar velocity 
dispersion we obtained weak correlation with the exclu- 
sion of low ionization lines, which are formed, presumably, 
in the outer regions of the NLR. Then it is reasonable 
that at large distances from th e nucleus the mot ion is gov- 
erned by the gravitational field l|Walsh et al.ll200"^ . We com- 
pared the FWHMi with the stellar velocity dispersion ob- 
tained by STARLIGHT. The Spearman rank correlation co- 
eflicients between FWHMstar = 2.35(t* and FWHMgas are 
listed in Table [71 with the usual meaning of the symbols. In 
S2 galaxies clear correlations exist except for [Nelll]A3869 
and HellA4686, while in S-I galaxies only the Balmer lines 
show a correlation with stellar kinematics, even if weaker. 
We note that in the S2 sample, the low ionization lines 
show a better correlation, except for [Ol]A6300 which is 
a rather weak line. In particular, we observed a stronger 
correlation for [Nii]A6584 (rs=0.69) th a.n for [Oni]A5007 
(ra=0 .5 3). This confirms the re sults bv iNelson fc Whittle! 
I 1996h . iBotte et all (|2005l ) and ICreene fc Hoi (|2005l ). that 
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[Oni]A5007 is not a good indicator of stellar velocity dis- 
pe rsion. A b etter indicator is [Nll]A6584, as already shown 
bvlHol l|2009l 'l. The diagrams of FWHMstar vs. FWHMgas are 
shown in Fig.[27]for both [Nll]A6584 and [Olll]A5007. 

In many cases, it was necessary to fit [Olll]A5007 with 
two components, a narrower (hereafter first component) and 
a broader component (hereafter second component). Fig. 1281 
shows the FWHM of [Olll]A5007 distributions for both the 
components. The values are systematically lower in the S2 
sample with respect to the S-I sample. The largest differ- 
ence is in the second component: the median value is 520 
km s~^ for S2 and 650 km s~^ for the S-I sample. 80 per 
cent of the sample has a FWHM of the [O ill] second com- 
ponent lower than 720 and 950 km s~^ respectively for S2 
and S-I. An accurate analysis of H/3 and [O ill] A4959, 5007 
profiles was performed. The asymmetry measure was ob - 
tained using the algorithm introduced by IWhittld (|l985h . 
Once the continuum level of the line is defined, the algo- 
rithm calculates the 10, 50 and 90 per cent of the area 
under the line profile and the corresponding wavelengths, 
from the blue side of the line. Then two line widths named 
respectively a=A5o — Aio and b=A9o — A50 are defined, so 
the asymmetry is calculated as A=(a-b)/(a+b). These mea- 
sures are based on the area, so they are less sensitive to 
noise or effects of instrumental resolution. The main issue 
is the definition of the line profile limits to fix the contin- 
uum level. A number of different baselines was used and the 
values of the resulting parameters were averaged. In order 
to define the line boundaries and the continuum we applied 
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Figure 27. The FWHM of [Oni]A5007 (up) and [Nn]A6584 
(bottom) are plotted against FWHM of the stellar component 
(FWHMstar). Symbols are Uke in Fig.[T2l 

the same method described in Section |3l Once evaluated the 
line boundaries, we considered three pixels at each side of 
the line (the pixel corresponding to the line boundary and 
the adjacent ones) in order to define 9 different baselines. 
A maximum limit of ±1500 km s~^ was imposed to pre- 
vent the non-convergence: in these cases the continuum value 
plus 2 or 3 times ac was taken into consideration, especially 
for [O lll]4959 and H/3, which have sometimes a relatively 
low S/N ratio. The distributions of the measured asymme- 
tries are presented in Fig. 1291 A clear indication of a positive 
asymmetry (blue wing) of the [Olll]A5007 line with respect 
to H/3 is found. Of course the H/3 asymmetry distribution for 
the S-I sample is totally different because of the broad com- 
ponent. However, the mean value for both samples is around 
A=0. In order to check the results, the asymmetry analy- 
sis was performed also for [O lll]4959. There is agreement 
between [Oin]4959 and [Oni]A5007 (Fig.[30l) although, the 
asymmetry measure is more difficult in case of [Oni]A4959, 
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Figure 28. The distributions of the FWHM of [Oin]A5007 for 
the first component (bottom panel) and for the second component 
(top panel). Colours and lines are like in Fig. 1131 



0.3 ^ . ^ . . . . ^ . . . . ^ . r 




0.5 



-0.5 0.5 

H/S asymmetry 

Figure 29. The distributions of the asymmetry for [Olll]A5007 
(top panel) and for H/3 (bottom panel). Colours and line styles 
are as in Fig. 1131 



because this line is weaker. The Spearman rank correlation 
coeflicients are r^^O.GO {t = 34.6, N=2131) for the S2 and 
rs=0.68 {t = 21.0, N=514) for the S-I sample. On the other 
hand there is no correlation between [Olll]A5007 and H/3. 
The [O lll]A5007 asymmetry parameter gives highly positive 
asymmetry for 68 per cent of the S2 and about 74 per cent 
of the S-I. The KS test gives a very little difference between 
the samples (D=0.08, A = 0.14, significance=l per cent). 
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Figure 31. An example of double components fitting; the ob- 
served spectrum (solid line), the narrow component (dashed line) 
and the broader component (dotted line). Residuals are plotted 
on the top of the diagram. 
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Figure 30. Asymmetry correlations, A5007 vs. A4g5g (bottom 
panel), and A5007 vs. Ajj/j (top panel). Symbols are like in Fig. ll2l 



This result supports the idea that when the asymmetries 
are present we are looking in the inner regions of the NLR 
where the kinematics is more complicated and chaotic, and 
probably the gas is not simply moving in the gravitational 
potent ial but it is driven by winds and decelerating out- 
flows (|WagneJll997l : iKomossa et al.ll2008l ). Finally, we did 
not find any correlation between these asymmetry param- 
eters and the luminosity, the FWHM of both components, 
and A(V^). 



5.5 Energy and Accretion rate 

With the ionized gas mass and velocity values, it is possi- 
ble to determine the kinetic energy due to turbulence, ther- 
mal and bulk motion. The third one would require velocity 
field, which are not available in our case. We could only 
use the asymmetric components (broader or second compo- 
nent) found in the [Oin]A5007 line to derive the order of 
magnitude of the kinetic energy due to bulk motion. These 
energies are de t ermin ed using the formalism introduced by 
iFu fc Stocktoiil (l2007l ). The turbulent energy can be esti- 
mated from the measured velocity dispersion cr(H/3), 
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-Etur = Ma 



5.0 X lO'^'^Afio (J50 erg 



(9) 



where (J50 — (j(H/3)/50 km s ^ and Mio the mass expressed 
in 10^" M0. The thermal energy is 



£th = -MkT/mp = 2.5 • lO^^Afio Ta erg 



(10) 



where T4 = T/ 10000 A' The bulk kinetic energy is given by 
Bbuik = \mV^ = 6.2 X lO^'Mio Vlo erg (11) 

where V250 is the bulk velocity in V^/250 km s~^ units. The 
bulk velocity is estimated from the difference in the peak 
velocities of the [O in]A5007 components (5p parameter), the 
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Figure 32. The distributions of turbulent energy (top panel), 
thermal energy (middle panel) and bulk energy (bottom panel) 
in erg 

. Colours and lines are as in Fig. llijI 

narrow component is assumed to have the systemic velocity 
(Fig.EIJ. 

Thermal and turbulent energy distributions are plot- 
ted in Fig. 1321 The distributions are similar for S2 and S-I 
with a logarithmic median value around 52.7 and 54.7 re- 
spectively for the thermal and turbulent energy. The ther- 
mal energy is negligible compared to the turbulent energy. 
In order to derive the bulk kinetic energy it is necessary 
to obtain the [O ill] A5007/H/3 ratio of the second compo- 
nent and then apply Eq.[S] to estimate the involved mass. 
Unfortunately, only in 14 spectra it was possible to fit two 
components in both [O in] and H/3. However, by plotting 
5p (H/3) vs. 5p ([Olll]), we found a good agreement between 
the fitted values, and the flux ratio of the second component 
is between 3 and 10 (Fig.l33l) , whic h is consistent with the 
as sumption of Fu fc Stockton! l|2006l ) and with the measures 
oflcid Fernandes et al. 1 (|200ll ). Then, by assuming 10 as a 
representative value of the [O lll]A5007/H/3 ratio, from the 
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Figure 33. The peak difference velocities of [O in] A5007 (km s"'^) 
are plotted against the same for H/3, for spectra with 2 compo- 
nents in both the lines (see text). 

luminosity of [Olll]A5007 second component and assuming 
low density A'o — 10 cm~^ we obtained a mass distribution 
with a median logarithmic value around 6.8 for both the 
samples. Hence the masses involved in the bulk motions are 
about 1/3 of the total mass. If the bulk motions are con- 
nected to outflows and we assume a typical dynamical time 
scale of 10^ yr (see below), the resulting rate of mass outflow 
is ~ 1 Mq yr~^. Finally, from the velocity peak difference, 
we derived the bulk kinetic energy by means of Eq.[Tl] The 
logarithm of median bulk kinetic energy values is around 
53.5. The distributions (Fig. 1321 bottom panel) show that 
the 80 per cent of the samples have values lower than 54.5, 
this means that the bulk kinetic energy is approximately 
lower than one order of magnitude compared to the turbu- 
lent energy. We stress that this could be an upper limit, 
since we assumed that all the flux of the second component 
is emitted by outflowing mass. If we assume that the outflow 
is directly driven by the black hole (BH), we can estimate 
the mass accretion rate necessary to sustain the outflow. The 
input momentum rate from radiation pressure is: 

jj = 1.3 X 10^^(r?/0.1)Afacc dyne (12) 

where t] is the radiative efficiency and Mace is the mass ac- 
cretio n rate of the BH in units of Mq, yr"'^ (|Fu fc StocktonI 
l200d ). The momentum rate of the moving bulk mass is given 
by: 

p ~ M X Vbuik/idyn dyne (13) 

where tdyn is the dynamical time scale. We estimated tdyn 
from the NLR radius 

tdyn ~ i^NLR/Vbulk S (14) 

obtaining the same distributions for both samples. The me- 
dian value is lO'^ yr, and 80 per cent of the samples have 
idyn < 3 X 10^ yr. Finally, by assuming 77 ~ 0.1 and making 
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Figure 34. Mass accretion rate distributions. Colours and lines 
are like in Fig. 1131 

Ea.ll2lequal to Eg . 1131 we obtained an estimate of the Mace 
distribution (Fig. l34|l . The median value is around 0.003 M© 
yr~^, 80 per cent of the samples have Mace < 0.03 Mq yr~^. 
These values are in agreement with the generally accepted 
accretion rates in Seyfert galaxies. 



6 SUMMARY AND CONCLUSIONS 

A sample of 5678 Seyfert galaxies was selected from SDSS- 
DR7 by applying the O123 diagram and their emission-lines 
were measured by means of a dedicated automatic code 
based on a nonlinear least square fitting method. Two com- 
ponents were applied to H/3 and Ha when necessary, and to 
[O lll]4959,5007 emission-lines in order to analyze their pro- 
files. The galaxies were divided into two groups on the basis 
of the Hq FWHM values and of the VO diagnostic diagrams, 
obtaining 2153 Seyfert 2 galaxies and 521 Intermediate- type 
Seyfert galaxies. Our results, obtained through a detailed 
spectroscopic analysis of the NLR of these two samples and 
a comparison between their physical properties, are in agree- 
ment with what is expected by the Unified Model and add 
new details to the general knowledge of the AGN properties. 

(i) The reddening for the gas, obtained by means of the 
Balmer decrement, is on average stronger in S2 (<A(F)> 
~ 1.5) than in S-I (<A(V)> ~ 1). According to the Uni- 
fied Model this is expected, because of the different line- 
of-sight orientation with respect to the torus axis. Indeed, 
iRodrfguez-Ardila. Pastoriza. fc Donzellil (|2000l ) found A(y) 
= 0.663 ± 0.345 for a sample of 16 Broad Lined Seyfert 1 
(BLSl). In addition, the reddening of the stellar component 
shows lower A(F) values, weakly correlated with those for 
the gas. This suggests that the extinction is mostly caused 
by dust associated to the gas of the NLR. 

(ii) Both the samples show good correlation between the 
lines having similar ionization potential. 

(in) The absence of [Ar lv]A471 1,4740 doublet and the 
values of temperature derived by [Olll], [On] and [Sll] line 
ratios, strongly support the idea that the NLR of Seyfert 
galaxies is a multi-phase ionized gas. In particular, it is likely 
made of a very low density medium (A^'c ~ 1 — 10 cm""^), 
where filaments or clouds with low density (Ac ~ 10^ cm"'') 
and high density (Ac > 4 x 10^ cm~'') are moving. The high 
density is expected to suppress [Ar iv] lines and at the same 
time to pump [Olll]A4363 by collisional processes. 

(iv) The [Fe vil] coronal lines are detected in ~ 4 per cent 
of the S2 sample, and ~ 19 per cent of the S-I sample. The 
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analysis of the [Fe vil]A6087/[O lll]A5007 emission-line ratio 
indicates higher values for the S-I sample and suggests that 
the coronal lines are formed d eep insi de the NLR. This result 
is in agreement with lNagao. T aniguc hi. fc Muravamal (I2OO0I ) 
who found log[Fe vii]A6087/[O iii]A5007 = -0.94±0.39, - 
1.33±0.43, -1.81±0.29 for BLSl, S-I and S2, respectively. 

(v) The ionization parameter U is in general significantly 
higher in S-I than in S2, showing average values of about 
—2.4 and —2.7, respectively. However, if we isolate the S2 
galaxies showing high ionization lines, as Hell A4686 and/or 
[Fe vil]A6087, the distributions of U for the two samples 
become similar. This result can be explained in case of a 
thinner torus which let the deeper and more strongly ion- 
ized NLR to be observe d. A comparison can be made with 
iFerland fc Netzerl l|l983l ) who estimated typical values of -2, 
-2.5 and -3.5, for SI, S2 and LINERs, respectively, by means 
of photoionization models. 

(vi) A good correlation between gaseous and stellar kine- 
matics is observed only in S2 galaxies, and only in case of low 
ionization and Balmer lines. This suggests that the kinemat- 
ics of the NLR is dominated by the gravitational potential 
of the galaxy only in the less ionized regions. Therefore, we 
confirm and point out that the FWHM of [Olll]A5007 line 
cannot be use as a surrogate of the stellar velocity dispersion 
in Seyfert galaxies. 

(vii) In addition, about 75 per cent of galaxies in 
both samples show asymmetries on the blue wing of the 
[Olll]A5007 line profile, suggesting the presence of radial 
motions, inflows and/or outflows, of the highly ionized gas. 
Also H/3 lines show asymmetries, which are nevertheless not 
correlated with those observed in [O ill] , with the exception 
of few cases (< 10 per cent). [O ill] A5007 profiles with strong 
blue wings have been also fo und in a sample of radio- loud 
quasars bv lBrothertoiil l| 19961 ). The author claimed that the 
H/? profiles show strong red wings. 

(viii) However, energy balance calculations, based on the 
estimated ionized gas mass, which is on average around 
3 X 10^ Mq for both samples, clearly indicate that the turbu- 
lent energy is the dominant component, with values in the 
range 10^^ — 10^^ erg. The kinetic bulk energy is an order of 
magnitude lower, while the thermal energy plays a negligible 
role. 
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